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Abstract 

The continuous wavelet transform using the Morlet basis may be symmetrized by 
using an appropriate normalization factor. The cone-of- influence is addressed through 
a renormalization of the wavelet based on power. The power spectral density may 
be deconvolved with the wavelet response matrix to produce an enhanced wavelet 
spectrum. The enhanced transform provides maximum resolution of the harmonic 
content of a signal. 

1 Introduction 

The continuous wavelet transform using the Morlet basis [H |2] has become quite popular for 
data analysis [3llll[5l[6l[71[8l[9l[l0l[IIl[T2l[13l[Tlj. There is some variety in the literature as to 
the assignment of the normalization factors, and we propose a rearrangement so as to produce 
a symmetric forward and inverse transform pair. The cone-of-influence is addressed through 
a renormalization of the wavelet amplitude which keeps its power constant for a given scale. 
The renormalized power spectral density may then be enhanced by deconvolution with the 
wavelet response matrix, yielding the maximum spectral resolution of the harmonic content. 
We conclude by discussing the utility of these algorithms and point out a recent application. 

The lack of a quantitative power spectral density has long hampered wider adoption of the 
continuous wavelet transform for data analysis. A mathematical engineer wants more than 
just a pretty picture — being able to give a numerical estimate to the power carried within a 
particular frequency band is of practical importance, and the use of a wavelet rather than 
Fourier transform allows that estimate to be time dependent. When the data is limited in 
duration, interesting features may be located in that region where wavelet truncation has 
become a significant effect. By renormalizing the power, the useful range may be extended 
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beyond the cone-of-influence, and perfect reconstruction holds for all but the edge-most 
locations. 

The breadth in scale of the wavelet response to a pure signal tone is a consequence of its 
localization in phase space, as the spectral and temporal resolutions are inversely related. 
Again, what a mathematical engineer wants is a precise identification of the frequency spec- 
trum, in which features are maximally resolved. By treating the continuous wavelet trans- 
form as a theoretical apparatus acting upon a signal, one may essentially calibrate the device 
given some basic assumptions on the form of the signal components. Minimizing the discrep- 
ancy between the continuous instant wavelet power and the convolution of the calibration 
matrix with the estimate then yields the sharpest resolution of the time-varying harmonic 
content of a signal. 



2 Normalization and central frequency 

We first discuss the symmetrization of the transform and its effect on the central frequency 
employed. One may write the usual Morlet wavelet [6l [13] at scale s = 1/ fs = In /u^ and 
offset t with unity sample rate At' = 1 using the parameter rj = (t' — t)/s as the product of 
a scale dependent constant C, a normalized window $, and a normalized wave G, 
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where Ui ~ 2tt is the central frequency of the mother wavelet at unity scale and zero offset, 
tpiQ^t') = TT^-'^/^e"* /2g«'^i* _ T]-^g window has an extent of — [sx to [sx and the wavelet a 
length of Nf = 2 where the parameter x = 6 defines the resolution width. The mother 

wavelet is normalized to unit energy J'^^\ip^Q{uj)\'^duj = 1 so that the Fourier transform 

q{suj). Pulling over from the denominator of the inverse 
transform a factor of the scale s and including a factor ^/2 which represents the transform 
response at negative scales gives a normalization Cg = \/2C^/s which produces a transform 
with some very desirable properties. 



(2) 



The mother wavelet now has a norm of 2, which we interpret as including the response at 
negative scales to negative frequencies, and that of a scaled wavelet is now 2/s^. With appeal 
to the photon, the energy of a localized wave is proportional to its frequency E^, oc u = s~^, 
thus its power (energy per time) should be proportional to its energy over its period, P,, oc 



The analysis of a (mean-subtracted) signal y(t) = i?e(ylfce*'^'=*) with duration Nt 
proceeds by the symmetric forward and inverse transform pair 



cwT(.,t)= J2^sAt')yit') 



ICWT(t) 



t' 
Re 



EECt(^')CWT(.,t')A5 



(3) 
(4) 
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Figure 1: CWT power spectral density (a), mean wavelet power (b), integrated instant power 
(c), and reconstruction (d) for a test signal with components of unit amplitude and periods 
5, 15, and 50 with rms power of 1.5. An unlabelled tick appears at the scale of the signal 
duration Nt = 300. Overlaying the PSD in (a) are the cone-of-influence (solid) and the 
cone-of-admissibility (dashed). 

for ip*{'r]) = V' ("■??), with perfect reconstruction within the cone-of-influence and quantitative 
agreement between the estimated power and the sum of the squared amplitudes of the 
signal components (c/. Equations (6) and (9) by Frick, et al [6j). Use of logarithmic scale 
spacing requires retention of the factor As. The root-mean-square power spectral density 
PSD(s,t) = |CWTp is normalized such that the integrated area of an isolated peak in the 
instant wavelet power IWP((s) = PSD(s,t) returns half the square of the amplitude Ak 
of the signal component, Pj-ms = The margins of the PSD determine the mean 

wavelet power MWP(s) = A^^"^ PSD(s, t) and the integrated instant power IlP(t) = 
PSD(s, t)As. The analysis of a test signal of duration Nt = 300 time units with signal 
components of unit amplitude and periods of 5, 15, and 50 is shown in Figure [H The cone- 
of-influence is marked with a solid line in (a), and the more restrictive cone-of-admissibility 
denoting the first wavelet truncation at a given scale is marked with a dashed line. A trough 
appears at the scale of the signal duration Nt in (b), beyond which we identify the extremely 
low frequency ELF region where s > Nt. Apparent is the loss of transform response in (c), 
where the IIP falls below the rms power Prms = 1-5, and the loss of reconstruction in (d) at 
the signal edge. 

The central frequency given by Torrence and Compo [13] to unify the Fourier period 
and wavelet scale, Ai/si = 47r/(u;i + a/2 + uf) = 1 yielding uji = 2n — l/47r, is no longer 
appropriate for our normalization. Using the same test signal, we consider transforms with 
central frequencies 2n — l/47r, 2tt, and 2it + l/47r and forward scalings of s~^, and 

appearing in the CWT. The top row in Figure [2] displays the instant wavelet power 
for a single central frequency at the center of the transform t = Nt/2, and the bottom row 
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Figure 2: Instant wavelet power (top row) and its gradient at the central peak (bottom 
row) for forward transform normalizations labelled by column and central frequencies of 
2tt — 1/Att (dash-dot), 27r (dashed), and 27r + 1 47r (solid). The test signal has components 
of unit amplitude and periods 5, 15, and 50. Each peak in (c) has an area of 0.5, which is 
equal to the rms power of the signal component. 

shows its gradient for all three central frequencies in the vicinity of the central signal peak; 
similar graphs obtain for the other peaks, noting that the forward scaling of in (b) and 
(e) corresponds to that recently proposed by Liu, et al [15]. Only the transform with scaling 
produces peaks with an integrated area equal to half the sum of squared amplitudes, 
and we note that the locations of its peaks coincide with the signal periods for the central 
frequency of uoi = 2tt + 1/Att. The response of the symmetrically normalized CWT is that of a 
theoretical apparatus whose point spread function preserves the area of a Dirac distribution 
representing the power carried by a pure signal component of infinite duration with constant 
amplitude and period. 

3 Renormalization 

We next introduce a renormalization which compensates for the reduction in response outside 
the cone-of-influence. The cone-of-influence is defined by the e-folding time, te = s\/2 
for the Morlet wavelet, indicating that region beyond which the response of the CWT is 
significantly affected by the wavelet truncation, which begins at the cone-of-admissibility. 
Various algorithms have been proposed for its rectification [HI El El [El [12] ; however, we have 
found that algorithms which alter the shape of the analyzing wavelet also affect its frequency 
response. Thus, we are led to proposing a simple renormalization such that for transform 
coefficients outside the cone-of-admissibility the wavelet is given a norm of 2/s^. For wavelets 
truncated by either edge of the signal, the window $^ is shifted by an offset r relative to 
an unshifted window $o defining the time span t'. The length of a truncated wavelet ipr is 
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Figure 3: RCWT power spectral density (a), mean wavelet power (b), integrated instant 
renormalized power (c), and reconstruction (d) for a test signal with components of unit 
amplitude and periods 5, 15, and 50. The MWP for the RAWT is shown in (b) as a dashed 
line. Reconstruction is maintained until the edge-most data points. 



Table 1: Ratio of mean integrated power to 



PSD 


CWT 


RCWT 


^max 


600 


300 


600 


300 


MIF/Prms 


0.95754 


0.95697 


1.002 


1.0004 



defined to be the lesser of the raw wavelet length or the signal length, A^,- = mm{Nt' , Nt) . 
The offset r(t) is determined from either the center of the signal or the location of the cone- 
of-admissibility, and the algorithm to keep everything aligned gets a bit complicated: for 
t' = max(0, [sx — and t'{k) = [—[sx, [sx] indexed by k from 1 to Nf, if r < then 

t' t'[l,mm{Nt>,Nt)] + r', else t' t'[max{l,Nt' - Nt + l),iVt'] - r'. The end result is 
simply to truncate either edge of the wavelet as necessary. Then for the amplitude of the 
truncated wavelet ipr = C^^t-G, with Cg ^ Cg^t = C*s -\/27s^|VvF define the renormalized 
continuous wavelet transform RCWT. 

Using the same test signal as above, in Figure [3] we display the analysis using the RCWT; 
the reconstruction in (d) is noticeably improved, and the power estimation in (c) is not as 
affected near the signal edge. The apparent increase in the IIP over the rms value represents 
we feel an aliasing in time, rather than scale, of the total power, as the mean discrepancy 
from the rms power ^^(IIP — Pims)/Nt = 5.6 x 10~^ is small. In Table [H we display the 
ratio of the mean instant power MIP = "^IZT ^^PSD/At to the rms power for the CWT 
and RCWT, considering also an integration over scale which stops at the signal duration 
Smax = Nt rather than Smax = '^Nf 
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Figure 4: RCWT power spectral density (a), mean wavelet power (b), integrated instant 
renormalized power (c), and reconstruction (d) for a test signal with components of time- 
varying periods around 5, 15, and 50 and squared amplitudes of .1, 1, and .5 respectively. 



Table 2: Comparison of FFT and RCWT rms power. 



Y^tyVNt 


EfFFT 


^^^^2 A.MWP 




0.79793 


0.79793 


0.80908 


0.80777 



As pointed out by Frick a/ [6], wavelet truncation also affects the admissibility condi- 
tion. One commonly subtracts from 0° = e«'^i»? the DC component ipi{^) oc e~'^i/^ = ~ 
10~^ of the mother wavelet so that the zero mean wave becomes B = Q^—dc- Then, the adap- 
tive wavelet transform is defined by 6 = 6° — dg^t) where d^ t = (0°) = 9°$,-/ is 
the weighted mean of the remaining wave 0°. Normalization as above with Cs,t then defines 
the renormalized adaptive wavelet transform RAWT of Reference |16j. In practice, we have 
found that the RCWT neglecting admissibility outperforms the RAWT, shown in Figure[3](b) 
as a dashed line, by a small but noticeable margin; the troughs between peaks are slightly 
deeper, and the reconstruction is slightly better. The reason, we feel, is that the adaptive 
admissibility condition alters the shape of the wavelet, hence its frequency response. 

The hallmark of wavelet analysis is its ability to track signal components with periods 
that vary in time. Considering now a test signal of the same duration, for periods 5, 15, and 
50 we adjust the squared amplitudes to be 0.1, 1, and 0.5 respectively (rms power of 0.8) 
and impose independent sinusoidal variation to the periods on the order of the duration. In 
Figure H] we display the RCWT analysis of such a signal, and in Figure [5] we compare the 
FFT power spectrum |tl7j with the MWP. The height of the peaks in the MWP appears to 
represent the relationship between the powers of the signal components better than the peaks 
in the FFT. Also, shown in Table [2] is the agreement between the mean squared amplitude 
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Figure 5: The FFT power spectrum (dotted) compared to the RCWT mean wavelet power 
(sohd). Locations of the discrete FFT values are marked (*). 

of the signal and the sum of the FFT power spectrum, which are quite close numerically to 
the value of the integrated MWP — the discrepancy decreases with increasing resolution in 
scale. For comparison, the integrated MWP for the CWT analysis of this signal is 0.77325, 
which is significantly less than the rms power. 

4 Enhancement 

With our transform now responding like a theoretical apparatus for measuring a signal's 
power spectral density, one may apply the techniques of resolution enhancement common 
in the analysis of experimental data |18]. First one writes the point spread function as the 
response matrix R{s, s') defined by the integrated power of a wavelet of scale s convoluted 
with a signal component of period A = s'. For this analysis we take the signal components 
to be Hann windowed cosine functions for the duration of the wavelet. 



which gave the best performance among the options. We note that here one is making an 
assumption on the form of the underlying signal elements whose composition represents the 
original signal. Then, for each IWP in the PSD, the enhanced instant power EIP is the 
solution to the equation 



found in a least-squares sense with non-negativity constraints. Note that it is the redundancy 
in scale of the CWT which provides the resolution enhancement of the EIP. In practice, we 
apply the enhancement to the peak PSD using R — > 2R. The effect is to replace broad 
peaks in the IWP with sharp spikes at the scale of the corresponding signal component, 
as shown in Figure [6] (a) for the IWP at the midpoint of the duration of the signal in the 
paragraph above. The reconstructed IWP (dashed) differs slightly from the original IWP 
(solid), sometimes exceeding it as no constraint has been placed on preserving the norm. 
In general, one's wavelet response may extend beyond one's region of calculation for signal 
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Figure 6: EIP (spikes) and IWP (solid) at the midpoint of the signal duration. The recon- 
structed IWP (dashed) slightly exceeds the original IWP. Panel (a) used a tolerance of 10~^, 
and panel (b) used a tolerance of 10"'^ to suppress the low frequency components. 

periods near either cutoff, and the enhancement procedure is capable of recapturing the lost 
(uncomputed) power. If one's application indicates the signal is bandwidth limited to that 
region well within the cone-of-influence yet far from the Nyquist scale, then enforcement of 
a norm-preserving constraint during the minimization is suggested. The enhanced power 
spectral density EPSD(s,t) is then defined simply as the collection of enhanced instant 
powers. Interestingly, one may effectively suppress the low frequency contribution by raising 
the tolerance of the minimization as in panel (b), where the EIP represents only the main 
signal components. 

In Figure [7] we show the EPSD of our most recent test signal, as well as the mean 
enhanced power MEP and integrated enhanced power lEP, which exceeds the previous IIP 
by a significant margin for which we have not yet found a correction factor. The variation in 
scale of the signal periods is well-resolved within the cone-of- influence, and a video scanning 
through the EIPs is available as an online supplement. Reconstruction from the EPSD is not 
yet well-defined; however, one may attempt a reconstruction using the phase of the RCWT 
and the original renormalized basis as shown in (d), which might not be perfect but does 
faithfully represent the original signal. For signal components within the cone-of-influence, 
the EPSD provides the maximum resolution in scale available from the RCWT. 



5 Conclusions 

The utility of these algorithms should be apparent to anyone familiar with one dimensional 
data analysis and power spectrum estimation. The extension of the renormalization prescrip- 
tion to mult i- dimensional wavelet analysis is straightforward; less so for the enhancement 
procedure. The symmetric normalization adopted here presents a power spectral density 
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Figure 7: Enhanced power spectral density (a), mean enhanced power (b), integrated en- 
hanced power (c), and reconstruction (d) for a test signal with components of time- varying 
periods around 5, 15, and 50. The increase in the lEP is attributed to the reconstructed 
IWP exceeding the original IWP. 

which behaves exactly as it should, where its margins may put on a physical basis, and 
provides perfect reconstruction without the introduction of an arbitrary factor. The mean 
wavelet power represents the relationship between signal components better than the fast 
Fourier transform, and the integrated instant power agrees with the rms power of the signal 
components. 

A recent application of algorithms very similar to the ones presented here is found in [SI 
[TB] . In that work, by addressing the power spectral density of the historical sunspot record, 
a relation is found between the level of solar magnetic activity and the temperature observed 
in central England. One also may consider its application to signal encoding, manipulation, 
and compression, providing an alternate basis for reconstruction. For temporally resolved 
power spectrum estimation, the symmetric wavelet periodogram has become quite a useful 
tool indeed. 

In summary, the continuous wavelet transform may be normalized to account for the 
response at negative scales, resulting in a symmetric forward and inverse transform pair 
with perfect reconstruction. It may then be renormalized to account for wavelet truncation 
by keeping a constant wavelet power across scale, where neglecting the admissibility condi- 
tion results in better performance, extending the useful range beyond the cone-of- influence. 
By deconvolution with the wavelet response matrix, the enhanced power spectral density 
provides the maximum resolution in scale of the harmonic content carried by a signal. 
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